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SUMMARY 


A flutter analysis employing the kernel function for three- 
dimensional, subsonic, compressible flow is applied to a flutter-tested 
tail surface which has an aspect ratio of 3-5> a taper ratio of 0 . 15 ., 
and a leading-edge sweep of 30°. Theoretical and experimental results 
are compared at Mach numbers from 0.75 to O. 98 . Good agreement between 
theoretical and experimental flutter dynamic pressures and frequencies 
is achieved at Mach numbers to 0.92. At Mach numbers from 0.92 to O. 98 , 
however, a second solution to the flutter determinant results in a 
spurious theoretical flutter boundary which is at a much lower dynamic 
pressure and at a much higher frequency than the experimental boundary. 


INTRODUCTION 


In order to predict successfully the flutter characteristics of a 
low-aspect-ratio, plate-like lifting surface, a flutter analysis must 
employ accurate, three-dimensional, aerodynamic forces. Because of the 
prominent role of flutter in modern airplane design, the development of 
practical procedures for calculating three-dimensional, unsteady, aero- 
dynamic forces is a problem of increasing importance. 

A promising approach to this problem involves direct consideration 
of the integral equation relating the lift and downwash distributions 
of oscillating finite wings. Basic to the solution of this equation is 
the evaluation of its kernel, a function which is essentially an aero- 
dynamic influence coefficient representing the downwash at some point 
on a lifting surface due to a unit aerodynamic load at any other point 
on the surface. In reference 1, the kernel function for oscillating 
finite wings in compressible subsonic flow was reduced to a form which 
could be conveniently evaluated. By using this form of the kernel, a 
subsonic lifting-surface method for calculating the forces on a har- 
monically oscillating wing of arbitrary plan form and deflection mode 
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was developed in reference 2. This method employed a systematic numeri- 
cal solution of the integral equation and was suitable for programing 
in high-speed computing machines. 

In reference 3 , the lifting-surface “method of reference 2 was 
employed in a modal- type flutter analysis. This flutter analysis appears 
to be well suited to the study of thin low-aspect-ratio surfaces at sub- 
sonic speeds, since it is capable of accounting for the effects of com- 
pressibility, finite span, and chordwise deformation. Few comparisons 
between flutter boundaries, calculated by this analysis, and experimental 
flutter boundaries are presently available. The purpose of the present 
paper is to present such a comparison. 

In the present investigation a tail surface having an aspect ratio 5 

of 3.5 and a taper ratio of 0.15 was flutter tested at Mach numbers 
from 0.75 to 1.06. This model had a leading-edge sweep of 30° and a 
3.5-percent-thick biconvex section. The natural vibration modes of the 
model involved significant chordwise deformation. Theoretical flutter 
boundaries were calculated, for Mach numbers from 0 to O.98, by using 
the analysis of reference 3 (with minor modifications). In the analysis, 
the flutter mode was approximated by a linear combination of the first 
three experimentally measured natural vibration modes of the model. 

Theoretical and experimental flutter boundaries are compared herein. 


SYMBOLS 


A 


A ij 


b 0 


b t 

f 

fi 

e 


aspect ratio 
panel aspect ratio 

generalized aerodynamic force — i- 

b 0 Z 

streamwise root semichord, ft 
streamwise tip semi chord, ft 
frequency of vibration, cps 
frequency of the 1th natural vibration mode, cps 

structural damping coefficient 
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h(x,y,t) 


m(x,y) 


x,y,z 


instantaneous deflection of point on tail surface in flutter 
mode, ft 

displacement of point x,y in the ith mode of vibration normal- 
ized with respect to the maximum displacement in that mode 
hi(x,y) 

imaginary part of complex unknown ft 
indices denoting a particular natural vibration mode 
reduced frequency, v/ v 
exposed semispan, ft 

dimensionless series expression for aerodynamic load distri- 
bution due to motion in jth vibration mode 

Mach number 

generalized mass associated with ith mode of vibration 
m(x, y)hjldS, slugs 
S 

local mass per unit area at point x,y, slugs/sc^ ft 
dynamic pressure, Ib/sq ft 

generalized force associated with ith mode of vibration, lb 
generalized coordinate in ith mode of vibration, q-je 1 ^, ft 

complex amplitude of generalized coordinate in ith mode, ft 

real part of complex unknown H 
surface area of tail, sq ft 
time, sec 

velocity of airstream, ft/sec 



0 


cartesian coordinates (see fig. k) 
angular chordwise variable (see fig. 5) 
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a LE 

A 


M- 




P 


leading-edge sweep angle , deg 

taper ratio 
panel taper ratio 

panel mass ratio, ratio of exposed panel mass to mass of a 
truncated cone of air having lower base diameter 2bQ, 

upper base diameter 2b^, and height Z 

generalized mass in ith mode of vibration referred to a den- 
sity parameter, 4itpZ 2 bQ (see eq. (ll)) 

density of airstream, slugs/cu ft 


a 


complex unknown of flutter determinant, 


I^Yd * ig) 

\ 03 / 


03 

CDi 


angular frequency, 2rtf, radians/sec 

angular frequency of ith vibration mode, 2itf^, radians/sec 
Dots over symbol indicate derivatives with respect to time. 


EXPERIMENTAL INVESTIGATION 


Model Geometry and Construction 

The experimental data, presented herein, were obtained from one of 
a series of all-movable, horizontal stabilizers which were flutter tested 
in the Langley 8- foot transonic pressure tunnel. The model used for the 
present investigation had an aspect ratio of 3-5, a taper ratio of 0 . 15 , 
a leading-edge sweep of 30°, and a 3 *5-P ercen "k-thick biconvex section. 
Structural and geometric details of the model are shown in table I and 
in figure 1. 

As shown in figure 1, the model was built up of 0.010-inch-thick 
duralumin laminations. These laminations were covered with balsa which 
was shaped to. the desired external contour. The outer surface of the 
model was covered with a thin plastic film. Incorporated In the model 
center section were two ball-bearing trunnions which located the hinge 
line at 40 percent of the mean aerodynamic chord. For the present 
investigation, however, the actuator rib (see fig. l) was clamped to 
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prevent its rotation about the hinge line. This clamping of the center 
section served to increase the amount of chordwise deformation involved 
in the natural vibration modes of the model. 


Physical Properties of Model 

The natural frequencies of vibration were determined by exciting 
the model with two electromagnetic shakers located fore and aft near 
the root chord. During these vibration tests, the model was mounted 
in the wind tunnel exactly as it would be mounted for the flutter tests 
which were to follow. Time-exposure photographs were employed to deter- 
mine the natural vibration modes of the model. These photographs were 
taken with a camera which was mounted outboard of the model tip and 
directed spanwise toward the root. The upper surface of one panel of 
the model was painted flat black and had thin, white, chordwise lines 
located at every one- tenth of the exposed semispan. With the model 
moving in one of its natural modes, a time-exposure photograph was 
taken. In the photograph, the chordwise white lines appeared as bands, 
the thickness of these bands being indicative of the amplitude of 
motion. The first three natural modes and their frequencies are pre- 
sented in figure 3 • 


Instrumentation 

The model was instrumented with electrical strain gages. Two 
groups of gages were used on each panel of the model. The first group, 
which consisted of bending and torsion gages bonded to both the upper 
and lower panel surface, was located near the 50-percent-chord line and 
the panel root. The second group consisted of torsion gages bonded to 
the upper panel surface and was located near the 50-percent-chord line 
at 70 percent of the exposed semispan. During the tests, a recording 
oscillograph was used to record the signals from the various gages. 
These records were used to determine flutter frequencies and the onset 
of flutter. The strain-gage signals were also fed into a cathode-ray 
oscilloscope in such a way that a Lissajous pattern indicated the start 
flutter. At each test point, the tunnel Mach number, stagnation 
temperature, and stagnation pressure were recorded by a punchcard 
readout system. 


Tunnel and Support System 

The - tests were conducted in the Langley 8-foot transonic pressure 
tunnel which is a single-return tunnel having a rectangular, slotted 
throat. In this tunnel, stagnation pressure and Mach number are 
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independently variable. Some details of the tunnel test section have 
been presented in reference 4. 

The model was mounted at an angle of attack of 0° on a long, 

5 . 5- inch-diameter cylindrical fuselage which extended into the subsonic 
region upstream of the test section. The fuselage was considered a 
rigid mount, since the fuselage mass was many times greater than the 
model mass. The measured fundamental vibration frequency of the support 
system was 4.3 cycles per second. 


L 

Flutter Tests 6 

1 

Experimental zero-lift flutter points were obtained at Mach numbers ' 

from 0.75 to 1.06. The procedure used in obtaining flutter at a given 
Mach number was to increase stagnation pressure gradually until flutter 
was obtained. After flutter was obtained, tunnel conditions were held 
constant momentarily and then Mach number and stagnation pressure were 
reduced as rapidly as possible in an effort to save the model from 
destruction. All flutter points reported herein were obtained from 
the same model. Experimental results are presented in table II. 


ANALYTICAL INVESTIGATION 


This section presents the method of flutter analysis used in the 
present investigation. Because of the similarity between the present 
analysis and that of reference 3, the method is described only in gen- 
eral terms. 

In the present analysis, the fuselage side was assumed to act as 
a reflection plane. Hence, the plan form which was analyzed is that 
plan form which results if, in figure 1, the center section is removed 
and the panels are joined along their respective root chords. This 
resultant plan form and the coordinate system used in the analysis are 
presented in figure 4 . 


Equations of Motion 

In the present analysis, the flutter mode is approximated by a 
linear combination of the model's first three natural (orthogonal) 
vibration modes; that is, 

h(x,y,t) = q^(t)h-^(x,y) + qgC^hg^y) + ^("^^(xjy) 


( 1 ) 
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where q i (t) = q i e ia) ' t is the generalized coordinate for the ith degree 
of freedom and h^(x,y) is the associated normalized natural mode shape. 

If the Lagrangian formulation is applied and the procedure of ref- 
erence 5 is followed, the resulting equation of motion for the ith 
degree of freedom is of the form 

Mi'di + CDi^qi = Qi (2) 


where Mi and are, respectively, the generalized mass and the 

generalized aerodynamic force for the ith degree of freedom. If the 
definition = q^ e io)t is employed, equation (2) may he written in 
the form 


1 



M i<li + Q ± = 0 
or 


(3) 


which corresponds to equation (6) of reference 3* 


Generalized Mass 


The generalized mass for the ith degree of freedom Is defined as 


Mi = 



( 4 ) 


where m(x,y) is the distribution of wing mass per unit area, and 

is the ith normalized, natural mode shape. In the present investiga- 
tion, both m(x,y) and hi were determined experimentally, as 

described in the section entitled "Physical Properties of Model." In 
order to evaluate Mi, the wing was divided into 27 mass elements, and 

the displacement of the center of mass of each element h^^) was 

taken from the appropriate mode shape. The generalized mass for the 
Ith mode was then computed as 



k=l 

where m( k ) is the mass of the kth element. 


( 5 ) 
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Generalized Aerodynamic Force 

As shown in reference 3, the generalized aerodynamic force may he 
expressed as 



bnpV^i 

-.2 2 
h n u> 



(6) 


In equation (6), Lj is a series expression for the aerodynamic load 

distribution due to motion in the jth natural mode. Basically, two 
steps are required in the calculation of Qq. The first step is the 

determination of L* and the second step is the evaluation of the 

J 


surface integrals 



Reference 2 describes a collocation 


procedure which is used to determine L y This procedure consists of 

representing the aerodynamic surface loading by a series expression that 
automatically satisfies the conditions of load at the wing edges. The 
downwash angles at various control points on the surface are then used 
as boundary conditions to determine the values of the arbitrary coeffi- 
cients in the assumed series. The locations of the control points used 
in the present investigation are shown in figure 5. After Lj is deter- 


mined, the evaluation of Qi is simply a matter of performing the surface 


integrations 



dS. 


In the present investigation, this inte- 


gration was performed numerically by using an 11-point Simpson's Rule 
in the spanwise direction and a Simpson's Rule in the chordwise 

direction. The numerical integration lattice is shown in figure 5* 


Flutter Determinant 

By substituting equation (6) into equation (3), the ith equation 
of motion may be written in the form 


1 



+ 


4rtpV^7 

Rj_bQ 2 o ) 2 



dS + 



= 0 


(7) 



VJi M On 


\J1 H ON f 


2W 


or alternatively, 


i m 




where 


+ p (^l A il + ^2 A i2 + % A i3) = 0 
*0^1 

(8) 

Ai J - b J, h * L J as 

S 

(9) 

ko = 

^ V 

(10) 

]_ WpZ^bo 

Pi " Mi 

(11) 


If equation (8) is written for the three degrees of freedom (nat- 
ural modes) under consideration, a set of three simultaneous homogeneous 
equations results. If these equations are to have a nontrivial solution, 
their coefficient determinant must vanish. 


Hence, 
k 2 


-(D 


Q, + 


'll 


A 


12 


l 13 


kQ^i 


kQ 2 M-i 


V 21 


kQ 2 ^2 


l 31 


- m • 

a 32 


v 22 


kQ 2 ^2 


ko 2 ^i 

a 23 

kQ 2 M2 


k 0 2 ^3 


k 0% 


1 _ g) 2 fi + 




= 0 ( 12 ) 


where the complex eigenvalue Q is defined hy 

,2 


a 


- (?) 


(13) 
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In this definition of n, o> is the unknown frequency and g is a 
damping coefficient that "becomes zero at the borderline condition 
between damped and undamped motion. 


Determination of Theoretical Flutter Boundaries 

In the present investigation, theoretical flutter characteristics 
were determined for Mach numbers M from 0 to 0.98 and mass ratios p 
from 20 to 100 by the following procedure. 

For each combination of M and p, the flutter determinant was 
solved for a range of reduced frequencies kQ. Since expansion of the 

flutter determinant results in a cubic equation In the complex unknown 
ft, each combination of M, p, and kQ yields three values of ft which 

must be examined to see whether they meet the requirements for a flutter 
point. By definition, 


ft = (1 + ig) = R + ii 

\ O) / 


therefore , 


0) 

ay 2 


— and 

/R 



At each M, flutter characteristics were determined from plots of g 
and u)/cjq 2 against in which p was a parameter. Flutter points 

were taken as being those points where g = 0 and > As 

trations, curves of g and against l/kg are presented in fig- 


ures 6 and 7, respectively. The resulting flutter characteristics were 
plotted as flutter boundaries of dynamic pressure q required for 
flutter and flutter frequency ratio as functions of M with p 


as a parameter. Theoretical flutter boundaries corresponding to wind- 
tunnel conditions were calculated for comparison with the experimental 
flutter points. 


DISCUSSION OF RESULTS 


The g against l/kQ and caja^ against 1/kQ curves of fig- 
ures 6 and 7, respectively, are presented to illustrate the behavior 
of the flutter determinant roots at various Mach numbers. In the 
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present investigation, the behavior of these roots was found to be simi- 
lar throughout the Mach number range 0 ^ M S 0.90. Therefore, curves 
are presented only for M = 0, 0.75, and 0.90. These curves show that, 
for the range 0 = M ^ 0.90, the only critical root is that one asso- 
ciated with the second natural mode. It should be noted, however, that, 
as the Mach number approaches 0.90, the third-mode root approaches the 
g = 0 line . 


At a Mach number of 0.92 the character of the determinant roots 
has changed. At this Mach number, the first- and third-mode roots are 
both critical. The first-mode root predicts flutter for the mass-ratio 
range from 20 to 100 and the third-mode root predicts flutter for mass 
ratios 60 to 100. From figure 6(d) it is seen that the third-mode root, 
loops up over the g = 0 line and then continues on the negative side. 
Hence, the third-mode root indicates flutter over a limited range of 
l/kQ and p, whereas the first-mode root indicates flutter for all 


V* 0 


greater than that value for which g = 0 and for 


20 5 p < 100. 


At Mach numbers of 0.94, 0.96, and 0.98, the first- and third-mode 
roots remain critical. The principal effect of an increase in M from 
0.92 to O .98 is that the flutter region predicted by the third-mode 
root grows until, at M = O. 98 , only a lower limit on 1/kQ is indi- 
cated within the reduced- frequency range investigated. 


Figure 8 presents theoretical flutter boundaries of dynamic pres- 
sure q required for flutter and flutter frequency ratio as 

functions of Mach number with mass ratio as a parameter. Figure 8 shows 
that the flutter points taken from the first and second-mode roots of 
the flutter determinant form smooth boundaries. These boundaries, which 
are shown as solid lines in figure 8, are determined by second-mode-root 
flutter points at 0 = M ^ 0.90 and by fir st-mode -root flutter points 
at 0.90 < M ^ O. 98 . As mentioned previously, the third-mode roots of 
the flutter determinant also yield flutter points at Mach numbers from 
0.92 to O. 98 . The flutter boundaries due to these points are shown, in 
figure 8, as dashed lines. In figures 6 and 7 it was seen that the 
third-mode roots of the flutter determinant indicated flutter over 
limited ranges of reduced frequency and mass ratio at Mach numbers 
from 0.92 to O. 96 . Hence, at these Mach numbers an upper and lower 
boundary of q and c ojo^ would exist for each p. In figure 8, only 

the lower boundaries are shown. The curves of figure 8 show that, 
throughout the range of variables investigated, flutter q for a given 
M varies with p (and hence with density). Therefore, these theoreti- 
cal results do not agree with the so-called constant-q concept of 
flutter. However, the variation of q with p becomes small at 
• high p. The effect of p on the variation of q with M is also 
of interest. Consider the first- and second-mode-root curves. For a p 
of 20, q decreases slowly as M goes from 0 to 0.90 and then increases 
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abruptly for M > 0.90. For p = 40, the decrease in q with increasing 
M is greater than for p = 20, and the minimum value of q is reached 
at M = 0.925. In general, figure 8 shows that, as p increases, the 
decrease in q with increasing M becomes more pronounced, and the 
minimum value of q occurs at a higher M. These trends bear a remark- 
able similarity to the trends of experimental data presented in figure 21 
of reference 5* 

In figure 9, the theoretical and experimental results of the present 
investigation are compared. The theoretical curves of figure 9 were 
obtained from calculations employing the experimental value of p for 
each M. Figure 9 shows that, for M from 0.75 to 0.92, the agreement 
between theory and experiment is good for both q and 0 aja^. For M 

from 0.92 to O. 98 , however, two theoretical flutter boundaries exist. 

In flutter prediction by theoretical means, the flutter solution which 
yields the lowest q must be accepted at each Mach number . Hence , on 
this basis, the third-mode-root boundary, which predicts a flutter fre- 
quency very near the frequency of the third natural vibration mode, 
must be accepted for 0.92 = M = O. 98 . For this boundary the agreement 
of both q and with experimental results is poor. If, however, 

the third-mode root were disregarded, the agreement between theory and 
experiment, especially with regard to trends, would be acceptable for 
the entire M range for which calculations were made. There is, how- 
ever, no reasonable basis for disregarding the third-mode -root flutter. 

The operating path of the wind tunnel was such that, if this flutter 
had existed physically, it would have been encountered during the « 

experimental investigation. 

In an attempt to determine the origin of the third-mode -root flutter, 
some calculations were undertaken for M = 0.94 in which various two- 
degree- of- freedom subcases of the flutter determinant were examined. 

When the first two natural vibration modes were considered, only the 
flutter solution which gave acceptable agreement with experiment was 
obtained. The third-mode-root flutter was found only when a combination 
of the second and third natural modes was used. It is noted that the 
frequencies of these modes are very close together j this in itself might 
be expected to lead to conditions of resonance. Since no such flutter 
as that indicated by the third-mode root was found experimentally, it 
may be concluded that this flutter mode is sensitive to some factor not 
properly accounted for in the analysis. A complete understanding of 
the nature of this additional flutter solution would require a broad 
survey of the many factors which might possibly influence it. Such a 
survey would include considerations of various control point locations 
and numerical integration techniques, of higher structural modes, and 
of structural damping. 
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CONCLUDING REMAP KS 


A flutter analysis employing the kernel function for three- 
dimensional, compressible flow was applied to a flutter-tested tail 
surface which had an aspect ratio of 3«5, a taper ratio of 0 . 15 , and 
a leading-edge sweep of 30°. The flutter tests covered the Mach number 
range from 0.75 to 1.06. Theoretical flutter boundaries were calculated 
for Mach numbers from 0 to O. 98 . In the calculations, the flutter mode 
was approximated by a linear combination of the first three natural 
vibration modes of the model. 

Good agreement between theoretical and experimental flutter dynamic 
pressures and frequencies was achieved for Mach numbers up to 0.92. At 
Mach numbers from 0.92 to 0.98, a second solution to the flutter deter- 
minant resulted in spurious theoretical flutter points which were at a 
much lower dynamic pressure and a much higher frequency than the experi- 
mental points. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., February 11, i 960 . 
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TABLE I 

PROPERTIES OF MODEL 


Aspect ratio, A 

Thickness ratio 

Taper ratio, A 

Leading-edge sweep angle, App, deg 

Streamwise airfoil section . . • 
Streamwise root semi chord, 1>q, ft 

Exposed semi span, Z, ft ■ . « » « 
Exposed panel mass, lh-sec^/ft 
Panel aspect ratio, Ap 

Panel taper ratio, Ap 

Natural frequencies, cps: 

First mode 

Second mode 

Third mode 


3.5 

0.035 

0.15 

30 

Biconvex 

0.625 

1.22 

0.055 

1.66 

0.176 

40 

108 

117 


TABLE II 

EXPERIMENTAL FLUTTER POINTS 


Point 

M 

V, 

ft/sec 

P > 

lb-sec^/ft^ 

<1, 

lh/sq ft 

f, 

cps 

O), 

rad/sec 

o)/oj2 


1 

0.750 

841.7 

0.002146 

760 

64.1 

402 

0.592 

44.3 

2 

.848 

938.5 

.001509 

664 

57 A 

360 

.531 

63.O 

5 

.912 

1,001.3 

.001197 

600 

51.9 

325 

.479 

79-5 

4 

.958 

1,044.1 

.001114 

607 

54.9 

344 

.507 

85.4 

5 

.983 

1,067.0 

.001232 

701 

56.5 

354 

.522 

77-2 

6 

1.060 

1,135.7 

.001114 

718 

56.1 

352 

.519 

85.4 
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Figure 1.- Structural and geometric details of model. All dimensions are in inches unless 

otherwise specified. 



Figure 2.~ Mass distribution of model 







Position of model 
At maximum amplitude 
At rest 



Figure 3 • - Experimentally determined natural mode shapes 




Position of model 



Figure 3 -- Continued. 



Position of model 
At maximum amplitude 
At rest 



Figure 3-- Concluded. 
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in inches. Note: z positive downward. 






Figure 6.- Continued 







First mode 
Second mode 
Third mode 



Figure 6.- Continued 









Figure 6.- Continued 






First mode 
Second mode 
Third mode 



Figure 6.- Continued 





Figure 6 . - Continued 






Figure 6.- Concluded 
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Figure 7*- Continued, 





Figure 7-~ Continued 
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Third mode 



Figure 7,- Continued 














